Confusing nonzero #13 with nonstandard interactions in the solar neutrino sector 
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Solar and KamLAND data are in slight tension when interpreted in the standard two-flavor oscil- 
lations framework and this may be alleviated allowing for a nonzero value of the mixing angle #13. 
Here we show that, likewise, nonstandard flavor-changing interactions (FCI), possibly intervening 
in the propagation of solar neutrinos, are equally able to alleviate this tension and therefore consti- 
tute a potential source of confusion in the determination of #13. By performing a full three- flavor 
analysis of solar and KamLAND data in presence of FCI we provide a quantitative description of 
the degeneracy existing between #13 and the vectorial coupling characterizing the nonstandard 
transitions between v e and v T in the forward scattering process with d-type quarks. We find that 
couplings with magnitude e^Y ~ 10%, compatible with the existing bounds, can mimic the nonzero 
values of #13 indicated by the latest analyses. 

PACS numbers: 14.60.Pq, 13.15.+g 



I. INTRODUCTION 

After many decades of efforts neutrino oscillations have 
been definitively identified as the leading mechanism gov- 
erning the flavor transitions observed in a variety of ex- 
perimental setups exploiting both natural and artificial 
sources of neutrinos. In analogy with the quark sector, 
where the mixing is described by the unitary Cabibbo- 
Kobayashi-Maskawa matrix, the leptonic mixing is also 
described by a nontrivial matrix connecting the flavor 
(y ai a = e,n,r) and the mass eigenstate neutrinos 
(i/j, i = 1,2, 3). However, it has long been noted [l[ that 
the matrix describing the propagation of light neutrinos 
is substantially more complex than the quark mixing ma- 
trix in that (i) deviations from unitarity may appear and 
(ii) there are additional CP phases, which have no quark 
analogue. Here we neglect both. Regarding the first we 
simply assume that their magnitude is small, as expected 
in high-scale seesaw schemes 1 . It is also well-known that 
Majorana CP phases affect only lepton-number violat- 
ing processes 0, [n| , and therefore can be neglected 
when discussing conventional neutrino oscillations. 

Hence, for simplicity, here we adopt the unitary ap- 
proximation for the lepton mixing matrix U describing 
neutrino oscillations, for which we take the standard fac- 
torized parametrization given in [l[ 

U = U(6 23 )U(9 13 )U(6 12 ). (1) 

as a product of three complex rotations characterized 
by three mixing angles (9 12, 6*13, 6*23) and three corre- 
sponding CP violating phases, and adopt the ordering 
prescription of the PDG [13]. There are two Majorana 



Low-scale seesaw models 

QHHHH would provide an exception 
to this, leading to large flavor and CP violating effects in the 
charged lepton sector [3, @, @, 0, fTojl , as well as effects in neutrino 
propagation lilt \lj Il3l . 



phases which do not contribute to neutrino oscillations; 
the Dirac phase is set to zero since current experiments 
show no sensitivity. 

The small mixing angle 613 is still unknown and its 
measurement constitutes one of the major goals in parti- 
cle physics, as it will open the door to possible measure- 
ments of CP violation in the leptonic sector [H, [l9| . 

Sensitivity studies of future experimental setups (20I 
l2lj | have evidenced how the identification of 6*13 can 
be problematic due to a potential confusion problem 
which may arise if nonstandard interactions (NSI) are 
present. These new interactions typically arise in low- 
scale models of neutrino mass, such as radiative ones [22I . 
I23I [2^ |. in the form of low-energy four-fermion oper- 
ators O a j3 ~ v a vpf '/, inducing either flavor-diagonal 
(q = 0) or flavor-changing (a 7^ 0) neutrino transitions 
in the forward scattering with the background / fcrmion 
dHHEHEHzl- More specifically, flavor-changing interac- 
tions (FCI) inducing transitions among v e and v T (a = e, 
(3 = t) have been recognized as an important source of 
confusion, as they can mimic the effect of nonzero #13 at 
neutrino factories 0, HH ■ 

The possibility that an analogous difficulty may be al- 
ready present in the interpretation of the available neu- 
trino data has not been considered so far. Indeed, in the 
past years the null result reported by the short-baseline 
CHOOZ reactor experiment [28[ (sin 2 6*13 < few%) was 
corroborated by the independent findings of the global 
neutrino data analyses. As a consequence, various works 
aimed at establishing the (subleading) role of NSI in the 
neutrino oscillation phenomenology have focused on the 
well-motivated two-flavor limit, in which #13 = is as- 
sumed. Furthermore, none of such analyses have evi- 
denced any significant preference for NSI. 

Recently, however, a nonzero value of 613 has been 
hinted in various analyses [H, M HE HI HI HI HS HH 
of the latest neutrino oscillation data. In particular, all 
the existing analyses [3l|, [H, Hf| find such a feature in 
the "solar sector" (solar and KamLAND data), while its 
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presence in the "atmospheric sector" (atmospheric and 
v n ~ * v [i disappearance long-baseline data) is more un- 
certain, being found in some analyses [H], [H, 
but not in others [33l . [38j . Interestingly, the preliminary 
searches performed by the Main Injector Neutrino Os- 
cillation Search (MINOS) in the v^ L — > v e appearance 
channel (39j seem to support such hints, showing a weak 
preference for a nonzero value of #13, just below the upper 
limit established by CHOOZ. 

In view of these hints we deem timely to investigate 
whether a confusion problem may already exist in the 
interpretation of the present data. Such an issue seems 
even more pressing, considering that the clearest hint 
of nonzero #13 comes from the solar sector, which natu- 
rally offers a sensitive setting, where the first signs of NSI 
may possibly emerge. Indeed, one should note that the 
"solar" hint of nonzero #13 arises from a tension among 
the standard interpretation of two-flavor transitions in 
matter (solar i/'s) — where NSI may intervene — and 
in vacuum (KamLAND) — where NSI are unimportant. 
Although the standard ?>v interpretation (#13 > 0) seems 
the most natural one, the possibility that this tension can 
be the result of some unknown effect intervening in solar 
flavor transitions cannot be discarded a priori. Here we 
investigate the possibility that such an effect may result 
from the theoretically well-motivated FCI, analyzing in 
detail their impact on the extraction of the estimates of 
6*i3 from the presently available data. 

The paper is organized as follows. In Sec. II, we review 
the notation necessary to describe three- flavor transitions 
in matter in the presence of NSI. In Sec. Ill, we present 
the results of our numerical analysis in the framework of 
2^ (6*i3 = 0) and Zv (#13 > 0) matter transitions in the 
presence of FCI. We draw our conclusions and discuss 
future perspectives in Sec. IV. 

II. NOTATION 

The Zv evolution in the flavor basis (v e , Vfi> v T ) is de- 
scribed by the equation 
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where H is the total Hamiltonian, 

H = H k[n + Hf^ n + iJ^yn 1 , (3) 

split as the sum of the kinetic term, the standard MSW 
(Mikheev-Smirnov-Wolfcnstcin) matter term [25|, [40j ] . 
and of a new, NSI-induced, matter term [ill ]. Indicat- 
ing with U the 3x3 mixing matrix, the kinetic term 
reads 

/ -5k/2 \ 
ffki„ = U\ +Sk/2 t/t , (4) 
\ k/2 J 



where E is the neutrino energy and 8k = Sm 2 /2E, 
k = m 2 /2E (8m 2 and m 2 being the "solar" and "atmo- 
spheric" neutrino squared mass differences, respectively). 
The second term Hf^ n describes the standard (MSW) dy- 
namics in matter [25l [4fj| , and is given by 

H?£ = diag(y, 0, 0) , (5) 

where V(x) = ^/2GfN c (x) is the effective potential in- 
duced by the interaction with the electrons with number 
density N e (x). The term characterizing the nonstan- 
dard dynamics, assuming for definitcness interactions 
only with d-type quarks, can be cast in the form 

(H^) aP = V2G F N d (x)e Q f3 , (6) 

where e a p = e^i are the dimensionless vectorial cou- 
plings between neutrinos with flavors (a, (3) with d-type 
quarks having number density Nd(x). In the phenomcno- 
logical approximation of one-mass-scale dominance, 

8m 2 < m 2 , (7) 

we can take the limit m 2 — > 00, and, similarly to the 
standard MSW case [4l[, reduce the 3v dynamics to an 
effective 2v one [42|. In fact, the 3x3 mixing matrix 
U(9i2, 6*13,6*23) can be factorized (assuming no CP vio- 
lating phase) into three real rotations 

U = R = R(6 23 )R(6 13 )R(6 12 ). (8) 

Performing a rotation of the initial neutrino (flavor) basis 
by R T (9i3)R T (623), and extracting the submatrix with 
indices (1, 2) one finds that the survival probability of 
solar electron neutrinos is given by 

p ee = 4?: ff + 4 , (9) 

where Sy = sin6>y , = cos6>y , and P^f is the v e sur- 
vival probability in an effective 2x2 model described by 
the Hamiltonian 

H cS = V(x) °) + V2G f N d (x) (° ^ , (10) 

where e and e' are two effective parameters which, con- 
sidering only FCI [a 7^ (3 in Eq. ©], are related to the 
original e a p couplings, as [i^] 

e = Ci 3 (e eAi C23 ~ £erS23) ~ SlS^T^ - S% z )] , (11) 

(' = -2e^ T C23S23 + 2si 3 Ci 3 (e er C23 + £e M S23) 

2 ( 12 ) 

- 2.s 13 e Mr s 2 3C23 ■ 

Considering we are interested in only the FCI between 
v e and v Tl we then remain with the expressions 

e = -e er ci3S23 i (13) 
e' = +2e er si 3 ci3C23 • (14) 

Therefore the propagation of solar neutrinos can be 
described effectively by a two-dimensional evolution 
Hamiltonian, which depends on the five parameters 

(Sm 2 1 e 12 ,e 13 ,e 23 , 
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FIG. 1: The solar LMA region is represented at two C.L.'s 
[Ax 2 = 1 (thin dashed line) and A% 2 = 4 (thick solid line)] 
for the standard case (ef^ = 0) and for two representative 
nonstandard cases with FCI having equal amplitude and op- 
posite sign (e^Y ~ ±0.2). The horizontally elongated regions 
are those allowed by KamLAND (at the same C.L's.) which 
do not depend on FCI. 




FIG. 2: Region allowed by the combination of solar and Kam- 
LAND data at two C.L's [Ax 2 = 1 (solid line) and Ax 2 = 4 
(dashed line)] after marginalization of 5m and #12. 



III. NUMERICAL RESULTS 

In our analysis we have included the data from the ra- 
diochemical experiments Homestake [43[ , SAGE [44[ and 
GALLEX/GNO [H, EE Super-KamioKande [H, 
from all the three phases of the Sudbury Neutrino Ob- 
servatory (SNO) Hi US El El, and Borexino jM- We 
have also included the latest KamLAND data [5J]. For 
the sake of precision, we have incorporated both standard 
and nonstandard matter effects in KamLAND. However, 
due to the low density of the Earth's crust, both have 
only a negligible effect for the range of parameters we 
are considering. Therefore, the constraints obtained on 
KamLAND do not depend on NSI. We also included NSI 
effects in the propa gati on of solar neutrinos in the Earth 
which, as noted in [55j |. can modify the regeneration ef- 
fect. 

We begin our study considering the more familiar two- 
flavor case (0i3 = 0) in which the results of our analy- 
sis depend on the three parameters: (Sm , $i2,£er 

). In- 
deed, we can safely assume sin 2 6*23 = 1/2, motivated 
by the results of the latest atmospheric neutrino data 
analyses, which indicate maximal [33[ or nearly maxi- 
mal (32l . [56j mixing. In Fig. 1 we show the region al- 
lowed by KamLAND (horizontally elongated region) in 
the plane spanned by the standard oscillation parameters 
[Sm 2 , sin #12] , superimposed to the solar large mixing an- 
gle (LMA) region of oscillation parameters obtained in 
the absence of FCI (e eT = 0), and for two representative 
cases in which FCI are "switched on," with couplings hav- 
ing the same amplitude but opposite sign (e er = ±0.2). 
The most important clue we have from this plot is the 
shift of the solar LMA region in the horizontal direction 



in correspondence of the (phcnomenologically relevant) 
values of 5m 2 determined by KamLAND. Indeed, the so- 
lar LMA region "moves" towards higher (lower) values of 
012, for positive (negative) values of e eT . From Fig. 1, it is 
clear that positive values of e eT tend to reduce the tension 
between solar and KamLAND and are preferred in the 
solar+KamLAND combination, which gives e eT ~ 0.15 
as best fit and disfavors the standard case (e eT = 0) at 
- 1.3cr level (Ax 2 ~ 1.7). 

These results strongly suggest that positive values of 
e eT , with size compatible with current limits (57j . can 
mimic the effect of nonzero #13. Notice, however, that 
in the standard analysis both the solar LMA re- 
gion and the region allowed by KamLAND in the plane 
[Sm 2 , sin 2 812] get modified b y va lues of 6*13 > 0, for which 
they tend to merge [la" [3fj|. l3ll. I33L |35| . In contrast, in 
the case under consideration, only the solar LMA region 
is affected. 

In order to trace more quantitative conclusions it is 
necessary to perform a full + FCI analysis, where 
both #13 and e eT are allowed to assume nonzero values. 
In Fig. 2 we display the main result of such an analy- 
sis, by showing the constraints (at Ax 2 = 1,4) obtained 
from the combination of solar and KamLAND, in the 
plane of the two relevant parameters [sin 2 6*13, e er ], af- 
ter marginalization over the remaining parameters Sm 2 
and sin 2 6»i2. We display only the region corresponding 
to positive values of e eT , which are relevant for the de- 
generacy problem under study. One sees that, at low 
confidence levels the preferred region is a band (delim- 
ited by the solid curves) which does not contain the origin 
("disfavored" at Ax 2 — 2.0). The slight tension among 
solar and KamLAND gets effectively diluted among the 
two parameters #13 and e eT . For higher confidence lev- 
els (dotted curve), the allowed region does contain the 



4 



origin and only upper limits can be put on both parame- 
ters. These results allow us to conclude that a complete 
degeneracy between the two parameters is present in the 
current neutrino data. 

The question arises as to whether and how future data 
may remove such a degeneracy. With this purpose, in 
Fig. 3 we show the behavior of the solar v e survival prob- 
ability (averaged over the 8 5 v production region) profile 
P ee (E), for three representative cases. In all of the three 
cases presented the probability is calculated for the fixed 
values of the leading parameters (5m 2 = 7.67 x 10 _5 eV 2 , 
sin 2 #i2 = 0.3). The three curves correspond to the fol- 
lowing cases: (I) The solid line represents the case of 
pure 1v standard transitions (#13 = 0, e eT = 0), corre- 
sponding to the origin in Fig. 2; (II) The dashed line 
indicates a representative case of standard ?>v transitions 
(sin 2 #13 = 0.02, e er = 0); (III) The dotted line shows 
a representative case of 2v + FCI transitions (#13 = 0, 

£ er =0.1). 

As one can infer from Fig. 2, cases (II) and (III) have 
been chosen so as to correspond to (currently) indistin- 
guishable points in parameter space. With respect to 
case (I), regarded as a benchmark, we note the following 
differences between the two degenerate cases, (II) and 
(III). In the standard 3v case (dashed line) P ee is sup- 
pressed with respect to the standard 2v case (solid line) 
by the energy independent factor ^1 — 2 sin 2 #13 [see 
Eq. JSJ], In contrast, the 2v + FCI case (dotted line), 
is characterized by an energy-dependent suppression 2 re- 
spect to the standard 2j/ case. In particular, the suppres- 
sion is completely negligible at low energies (E < 3 McV), 
and it is more pronounced at intermediate energies. The 
net effect is a flattening of P ee (E) with an enhancement 
of the up-turn typical of the adiabatic MSW transitions. 

The current data are unable to distinguish between 
case (II) and (III) since: (i) the differences at low energies 
are too small to be detected by the gallium experiments 
or Borexino; (ii) the current high energy experiments are 
not sensitive enough to probe the up-turn region, which 
still remains practically "invisible." The differences at 
low energy between the two degenerate cases, (II) and 
(III), are very tiny and may prove very hard to detect 
even at future low-energy experiments. Instead, the pos- 
sibility to disentangle the different behavior at interme- 
diate energies could perhaps become realistic in high- 
energy experiments with a lowered threshold. In this 
respect, the new data expected from Borexino, Super- 
K-III [531, and from the low energy threshold analysis 
underway in the SNO collaboration [59}, may play an 
important role. 

We close this section with a final remark. For definite- 
ness, we have focused on the case of interactions with d- 
type quarks and transitions among v e and v T . However, 
the essence of our conclusions is unaltered if interactions 
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FIG. 3: Solar v £ survival probability (averaged over the 8 B v 
production region) for three representative cases. 

with u-type quarks or electrons and/or transitions among 
v e and (a = e, /3 = fi) are considered. Moreover, the 
simultaneous inclusion of more than one type of NSI can 
only exacerbate the confusion problem we have posed, 
leading to further difficulties in data interpretation. 

IV. CONCLUSIONS 

We have stressed that nonstandard flavor-changing in- 
teractions may constitute a source of confusion in the 
interpretation of present solar and KamLAND neutrino 
data, hindering the correct determination of the mixing 
angle #13. In the near future, various solar experiments 
may help to reduce such a difficulty by providing a more 
precise determination of the energy profile of the solar v e 
survival probability. In the mean time, it would be very 
important to complement our study — focused on the 
solar sector — with a similar investigation on the atmo- 
spheric sector, which could hopefully be of aid in breaking 
the degeneracy among #13 and e eT . Also, a quantitative 
assessment of the impact of nonstandard neutrino inter- 
actions in the interpretation of the preliminary MINOS 
data in the — > v e appearance channel (39l | would be 
highly desirable. Our results underline the importance 
of a "clean" measurement of #13 expected from the new 
generation reactor experiments [601 ] , whose inferences are 
free from NSI effects. We also stress how, in the event 
of a null result by these experiments (i.e., in the case 
of nonconfirmation of the present hints of nonzero #13), 
a persisting tension among solar and KamLAND would 
pose a novel problem, whose resolution may involve the 
nonstandard interactions discussed here or other possible 
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unaccounted effects. 
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